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The equilibration of the right-helicity states v+ of light Dirac- neutrinos 
is discussed. I point out that the v + production rate is enhanced by 
■ weak gauge boson pole effects so that the right-helicity component 

Q\ . of v T is brought into equilibrium at T ~ 10 GcV independently of 

the initial abundance, provided m v > 10 keV. Neutrino spin flip 
in primordial magnetic fields and the resulting bound on fj, u is also 

D ' 

jy^ ' discussed. 

in 
cm . 

1 Relic abundances of right-helicity neutrinos 

Primordial nucleosynthesis is a remarkable probe of neutrino properties 0. To some 
extent primordial nucleosynthesis could be sensitive even to the Dirac vs. Majorana 



> 

nature of neutrinos, because in the Dirac case the small relic abundance of the inert 
right-handed component of Dirac neutrino would also contribute at nucleosynthesis. 
Of course, presently one cannot hope to differentiate between the Dirac and Majorana 
nature of neutrinos on cosmological grounds, but in principle this is an interesting 
Q_i! problem. The right-helicity states of Dirac neutrinos can be produced (and destroyed) 



in spin-flip transitions induced by the Dirac mass!'! or the neutrino magnetic moment. 

Spin-flip transitions may also be induced by primordial magnetic fields, if such exist. 
• i-H , 

^ ■ The actual cosmological density of the right-helicity neutrinos depends not only 

on the production rate near the QCD phase transition, but also on whether the right- 
helicity neutrinos had a chance to equilibrate at some point during the course of the 
evolution of the universe. This depends on the thermal scattering rates of neutri- 
nos. An important source of iz+'s are also the non-equilibrium neutrino scatterings 
and decays of pions!. Such processes give rise to the bound! m v < 130 keV and 
m VT < 150 keV, using Tqcd = 100 MeV and assuming that nucleosynthesis allows less 
than 0.3 extra neutrino families. Let us note that there seems to be no window of 
opportunity! for a sufficiently stable {r v > 10 2 sec) tau neutrino in the MeV region 
because of the production of non-equilibrium electron neutrinos in v T v T annihilations, 
which would disrupt the succesful nucleosynthesis predictions. 

In the Standard Model there arei 68 purely fermionic 2^2 processes in which 
a right-helicity muon or tau neutrino can be produced. In addition, a right-helicity 
tau neutrino can also be produced in 11 lepton and quark three-body decays, and the 
muon neutrino in another set of 11 three-body decays. In principle, one has also to 
consider processes involving W ± , Z and H. There are 16 such processes. Finally, 
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Figure 1: Thermally averaged production rate for the s-channel reaction ud — > v + t + with m VT = 20 
keV; shown is also the corresponding t-channel rate as well as the rate for the bosonic process 
t~7 — » v + W~ . 

there are 3 two-body decays of W ± , Z and H bosons which are capable producing 
right-helicity muon and tau neutrinos. 

The thermally averaged production rate in 2 — > 2 scattering a + b — > u + + d per 
one right-helicity neutrino v + is 

T+ = 4d f dU a dU b dU + dU d (27i) 4 5^(p a + p b -p + - Pd )S\M ab ^ +d \ 2 

77^_|_ J 

xr/ni-Hi-/? 5 ), a) 

where n+ D is the the equilibrium number density of the right-handed neutrinos, dUi = 
d 3 pi/((27r) 3 2Ei), S is the symmetry factor taking into account identical particles in 
the initial and/or final states, and ff B are Fermi-Dirac distribution functions. The 
fermionic processes exhibit an enhancement of the u + production rate which is due to 
gauge boson pole effects i. This is demonstrated in Fig. 1, where the enhancement 
is apparent in the s-channel process ud — > v + t + , where v + is a right helicity tau 
neutrino, as compared with the crossed t-channel process. It turns out that in the 
temperature of interest, the processes involving gauge or Higgs bosons can be neglected 
in comparision with the purely fermionic processes. This is demonstrated in Fig. 1, 
where the thermally averaged rate for the bosonic process r~7 — > v + W~ is shown. 

Let us note that at high T the rate Eq. ([!]) is infrared sensitive to the thermal 
corrections in the propagators. In most cases it is an excellent approximation just to 
modify the propagators by introducing a Debye mass M 2 (T) = U.L(uj,k = 0), which 
may be approximated by Mf (T) ~ Mf + 0.1 T 2 (i = W, Z). 

The relic density of the right-helicity tau neutrinos can be found by solving the 
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Boltzmann equation 




H \p + \»un)f + =\-£) ■ ( 2 ) 

coll 



Considering only purely fermionic 2^2 and 1 — > 3 processes, we can write the 
collision term in the form 

fdfA 



dR 



| = (C 2 ^ 2 + CU,)(1 - /+) - (C£_* + CU)f + , (3) 

coll 

where the coefficients Cj represent production and C' T destruction of u + , with I = 2 — ► 
2, f — > 3. The explicit expressions for the quantities C 2 ^ 2 and Ci^ 3 are 

C2^ 2 (|p + |,i2) = 2 tt4- /rfn a rfn feC /n d (2 7 r) 4 ( 5( 4 )(p a + j9 b -p + -p d ) 

„„„« /£/_!_ J 



scatt 



x^ QM i 2 /r/r(i-/r) 



dec + 

xS\M f „ g+h \ 2 f™(l-f™)(l-fn- (4) 

The Boltzmann equation Eq. (0) can be solved numerically (by using e.g. 30 
momentum binsi). One can consider two extreme initial conditions for v + at T ~ 100 
GeV: (i) full equilibrium with / + = / eq and (it) complete decoupling with / + = 
0. The evolution of the v + energy density p + for different masses and for the two 
different initial conditions is shown in Fig. 2. One observes that right-helicity tau 
neutrinos equilibrate independently of the initial condition provided m Vr > 10 keV. 
The relic density, in units of extra neutrino species AN„, is also shown in Fig. 2. 
For ra Vr < 10 keV the relic density is very close to the expected value p+/pi = 
[g*(T)/ g^Ti)} 4 / 3 p + (Ti)/ p L (Ti), where pl is the density of left-handed neutrinos and Tj 
is the initial temperature. For itl Vt > 10 keV one would obtain a small extra effect 
at nucleosynthesis, independently of the initial abundance. The results are for u T , but 
for the graphs would look very similar. 

2 Effects of magnetic fields 

The abundance of the right-helicity states is also sensitive to primordial magnetic 
fields. Large magnetic fields could be created by cosmological phase transitions at 
microscopic length scales i'l. There are indications that the small scale field cascades 
into larger length scales because of hydromagnetic turbulence!. Electrical conductivity 
in the early universe is hif?h0 so that the field is frozen in the plasma. Typically, such 
a primordial magnetic field is random so that (B) = 0. When a Dirac neutrino 
propagates through such magnetic field, its spin will precess. At the same time, the 
neutrino is subject to scattering and thermal corrections. Neutrino dispersion relations 
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Figure 2: (a) The evolution of the right-hclicity tau neutrino relic density p + , in units of the equi- 
librium density p FD ; (b) the relic density of the right-hclicity tau neutrino in units of extra neutrino 
species; the upper set of curves is for right-hclicity neutrinos in equilibrium at T ~ 100 GeV, the 
lower set corresponds to zero initial density. 
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ignetized plasma have recently been carefully studied^. The actual spin evoluti 
is best described by a relativistic kinetic equationlll, and the result depends very much 
on whether the coherence length of B is larger or smaller than the neutrino scattering 
length. Another unknown is how to average over the randomly varying magnetic field. 
In the case of neutrino propagation, the appropriate statistical procedure might be the 
line average!. Purely phenomenologically, one can write 



£ rm s = 5o(£) (£ 



(5) 



where Lq is the coherence length of the magnetic field, d is the distance scale and p is 
unknown; for a line average, p — 1/2. 

In the case of small-scale random magnetic field, neutrino forward scattering tends 
to depolarize the spin El. If the field is completely uncorrelated at distances d ^> Lq 
so that 

(B(t)B(t 1 ))=B* ms L 5(t-t 1 ) , (6) 
one finds that the spin flip probability reads 



1 



l-exp(-rt)) 



(7) 



where the damping parameter V is given by 



v — -n 2 R 2 T 2 

1 ~ gFv^rms-H) 



Setting t = H~ l and requiring that r < H at T = Tqcd so that the right-helicity 
states are not in equilibrium below QCD phase transition, and in particular during 
nucleosynthesis, results in the bound 

^B^Tqc^H- 1 ) < 6.7 x KTVbG f^V /2 • (9) 



In order to translate this to a bound on \x v one needs to assume something about the 
magnitude of B rms at the horizon scale at T = Tqcd- No reliable estimate exists at 
present time. However, the bound Eq. (|9]) could in principle be as restrictive as 
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